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We improve and extend Shapiro's Jl] model of a relativistic, compact object which is stable in 
isolation but is driven dynamically unstable by the tidal field of a binary companion. Our compact 
object consists of a dense swarm of test particles moving in randomly-oriented, initially circular, 
relativistic orbits about a nonrotating black hole. The binary companion is a distant, slowly in- 
spiraling point mass. The tidal field of the companion is treated as a small perturbation on the 
background Schwarzschild geometry near the hole; the resulting metric is determined by solving the 
perturbation equations of Regge and Wheeler and Zerilli in the quasi-static limit. The perturbed 
spacetime supports Bekenstein's |l3| conjecture that the horizon area of a near-equilibrium black 
hole is an adiabatic invariant. We follow the evolution of the system and confirm that gravitational 
collapse can be induced in a compact collisionless cluster by the tidal field of a binary companion. 
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I. INTRODUCTION 

The possibility that massive neutron stars might be 
driven unstable to collapse to black holes when placed 
in a close binary orbit was first suggested by Wilson, 
Mathews and Maronetti (hereafter WMM; [§) on the 
basis of their approximate relativistic numerical simu- 
lations. This finding was quite unexpected, partly be- 
cause it disagreed with earlier Newtonian calculations 
H which showed that tidal fields stabilize binary stars 
against radial collapse. In fact, none of the follow-up 
post-Newtonian (PN) |Q || or approximate analytic 
analyses || indicate the presence of any relativistic ra- 
dial instability in fluid binaries, nor does an indepen- 
dent dynamical simulation M. No evidence is found for 
the WMM "crushing effect" in the relativistic numeri- 
cal calculations of either corotational J^] or irrotational 
|)| binaries in quasi-equilibrium circular orbits. These 
later calculations are particularly careful to adopt the 
same simplifications as WMM (e.g, a conformally flat 3- 
metric). All of the calculations except those of WMM 
suggest that the maximum allowed rest mass of a fluid 
star in a binary is in fact slightly larger than the value 
in isolation. 

While it would appear that the WMM effect may not 
be present for fluid stars, the work of WMM raises the 
interesting question as to whether the "crushing insta- 
bility" might exist for a different type of binary sys- 
tem. Specifically, are there any relativistic binary sys- 
tems for which tidal fields can trigger the collapse of a 
compact object known to be stable in isolation? To ad- 
dress this question, Shapiro JlJ offered a simple candi- 
date compact object consisting of a test particle orbiting 
a Schwarzschild black hole outside the innermost-stable 
circular orbit (ISCO). Such a "compact object" is obvi- 
ously stable in isolation, but when it is placed in a binary 
orbit about a distant mass, an integration of the test- 
particle equations of motion reveals that the tidal field 



of the distant mass can cause a test particle to plunge 
inside the hole. 

Shapiro's model is very promising as a simple illustra- 
tion of the crushing instability in action, but the origi- 
nal analysis, as he emphasized, was highly simplified and 
somewhat heuristic. First, the study was confined to 
two special test-particle orbits, one coplanar with and 
the other perpendicular to the companion orbital plane. 
Moreover, the motion of the particle in the perpendicu- 
lar case was artificially constrained to remain in a plane, 
precluding any precessional motion or wobbling. Most 
important, the analysis was based on a post-Newtonian 
treatment of the 3-body problem in which the tidal piece 
of the equation for the relative motion of the test parti- 
cle about the black hole was treated to lowest (Newto- 
nian) order and the nontidal piece replaced by the ex- 
act, fully relativistic expression for geodesic motion in 
Schwarzschild geometry. While such a hybrid approach, 
similar in spirit to one proposed by Kidder, Will and 
Wiseman Q for the 2-body problem, should capture 
the essential dynamics, it is far from rigorous. In par- 
ticular, it is not clear a priori how reliable it is to use a 
Newtonian tidal term when the test particle is moving at 
high velocity in a very strong gravitational field close to 
a black hole. 

In this paper we improve and extend the simple model 
presented in [0 of a compact object subject to binary- 
induced collapse. We use the equations of Regge and 
Wheeler and Zerilli jlj] to derive the perturbation 
to the spacetime close to a Schwarzschild black hole due 
to a distant mass. We note in passing that the resulting 
perturbed spacetime provides another example support- 
ing the recent conjecture of Bekenstein 113] that the 
horizon area of a perturbed black hole is an adiabatic 
invariant. We obtain the geodesic equations of motion 
in the perturbed spacetime and solve them to track the 
dynamical evolution of a dense spherical swarm of 20,000 
test particles placed around the hole outside the ISCO. 
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The 'black hole + swarm' constitutes a compact colli- 
sionless cluster whose stability against tidally-induced 
collapse we determine. We compare our refined treat- 
ment with the simpler version presented in jjj and show 
that the original equations track the behavior reasonably 
well. We confirm that the crushing instability can occur 
in a binary system containing simple compact collision- 
less clusters like the ones we construct. 



II. HYBRID - PN TREATMENT 

A. Basic Equations 

Consider first a system of three bodies. Assume that 
two of the bodies, 1 and 2, say, are much closer to each 
other than they are to the third, so that the influence of 
body 3 on the relative orbit of 1 and 2 may be treated 
as a small tidal perturbation. Specialize to the case in 
which body 1 is a Schwarzschild black hole of mass M 
and body 2 is a test particle (nit "C M). Since it is far- 
away from the 1-2 pair (the "compact object"), body 3 
can be treated as a point particle, with mass m. Let r 
be the coordinate position of the test particle relative to 
the hole and R be the coordinate position of the distant 
point mass relative to the hole. Following Shapiro 0, 
we may write the equation for the relative motion of the 
test particle about the hole as 



r = — 



[Ae r + Bv} + — 



3R r 
~R?~ 



R r 



(1) 



where e r = J^. 

In the Newtonian limit, A = 1 and B = in Eq. (Q). 
For isolated binaries, Lincoln and Will Q] derive post- 
Newtonian expressions for A and B for arbitrary masses, 
correct through 2.5PN order. Kidder, Will and Wise- 
man fic| provide a "hybrid" set of equations in which 
the sum of the terms in A and B that are independent 
of the ratio n = /J./(M + nit), H = Mnit/(M + m f ), is re- 
placed by the exact expression for geodesic motion in the 
Schwarzschild geometry around a body of mass M, while 
the terms dependent on n are left unaffected. Their re- 
sulting equation of motion is therefore exact in the test- 
body limit (n — > 0) and is valid to 2.5PN order when 
appropriately expanded for arbitrary masses. We shall 
utilize these same hybrid expressions for A and B and, 
for simplicity, work in the test-body limit by taking one 
member of our close pair to have a mass much smaller 
than the other. Adopting harmonic (de Donder) coor- 
dinates, the resulting (Schwarzschild) expressions for A 
and B are given by 



1 - M/r 




2- M/r 


(1 + M/r) 3 _ 




1 - (M/r) 2 _ 



M 



•2 , 2 

r + v , 



B = — 



4 - 2M/r 



1 - (M/rf 



where 
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To model the three-body system, we must also know 
the position of the companion (body 3) relative to the 
compact object. The leading Newtonian piece of the 
equation of motion for R gives 



R 



M 



R 3 



-R. 
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We are interested in the dynamical behavior of the 
close pair, regarded as a single "compact object", as it 
inspirals toward the distant mass m due to gravitational 
radiation emission. To treat the inspiral of this "binary" 
(m in orbit about the "compact object"), we must include 
radiation reaction terms in the lowest order (Newtonian) 
orbit Eq. (^ . Formally, such a treatment requires a con- 
sistent expansion up to 2.5PN order. In lieu of this, we 
shall analyze the inspiral by assuming that the binary 
is in a nearly circular, Keplerian orbit, which undergoes 
a slow inspiral due to gravitational radiation loss in the 
quadrupole limit. This assumption is equivalent to in- 
serting a quadrupole radiation reaction potential in the 
binary orbit Eq. (j5j) and neglecting the lower-order, (non- 
dissipative) PN corrections in that particular equation. 
While such an expression is not formally consistent to 
2.5PN order, it faithfully tracks the secular inspiral of 
the binary in the limit treated here in which the binary 
system is at wide (nonrelativistic) separation. The de- 
tails of the inspiral are not important here, only that the 
inspiral serves to bring a tidal perturber slowly in from 
infinity toward our compact object ]l5| . The resulting 
equations for the binary inspiral are then JlJ 



R(t)/R(0) = (1-t/T) 1 ^ 



and 



e(t)-e(o) 



8 (M 



5 V R(0) 3 



1/2 



T[l- (1-t/T) 



5/8] 
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(Note the typo in Q.) The binary inspiral timescale T 
is given by 



T/M = 

We will use Eqs. 
Eq. (Eh. 



(R(0)/M) 4 



256 ((m + M)/M)(m/M) ' ^ 
(H)-(|l) in analyzing the relative orbit 



(2) 



B. Numerical Implementation 

To assess the fate of the compact object, we followed 
the motion of a spherical swarm of 20,000 test parti- 
cles about the black hole. At t = 0, the test parti- 
cles were placed randomly about the hole at a radius 
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r/M — 5.9, well outside the ISCO of an isolated black 
hole at r/M = 5, and set in circular orbits of arbitrary 
orientation. We set m = M and started the simulation 
with the companion at R/M — 60. At this initial sep- 
aration, the tidal field of the companion is negligible at 
the compact object, and the test particles begin in sta- 
ble equilibrium orbits. We integrate ([!]) until the time 
at which R/M = 26.30. Cartesian coordinates were used 
throughout. 

The results of the simulation are summarized in Fig- 
ure 1, where we plot the mean cluster radius as a func- 
tion of companion separation. Clearly, the compact ob- 
ject remains stable until the companion comes within 
R/M ps 30, at which point the tidal field causes many 
of the particles to plunge into the hole. We find that 
17,919 particles (89.6 %) have fallen into the black hole 
by t/M = 121,893, the time at which the companion 
reaches R/M = 26.30. The simulation confirms the exis- 
tence of a crushing instability. 

Several checks were performed to ensure the accuracy 
of our code, which integrates the ordinary differential or- 
bit equations by a standard fourth-order Runge Kutta 
algorithm with an adaptive stepsize ]l6[] . For example, 
in the absence of the companion, the code conserves par- 
ticle energy and angular momentum over the full time pe- 
riod of integration and reliably locates the Schwarzschild 
ISCO with a controllable precision. These tests are not so 
trivial when the integrations are performed in Cartesian 
coordinates. 



III. SCHWARZSCHILD PERTURBATION 
TREATMENT 

A. Basic Equations 

We will now improve the previous treatment by incor- 
porating the tidal field of the companion star in a fully 
relativistic fashion. To do this, we treat the effect of the 
companion as a small perturbation on a Schwarzschild 
background. By assuming that the companion orbits at 
a much greater distance from the black hole than the test 
particles, and hence moves with a much smaller angular 
velocity, we can work in the quasi-static approximation. 
We place the companion on the z axis of our spherical 
polar coordinate system centered on the black hole (see 
Figure 2). We follow Regge and Wheeler (ll) and divide 
the linear metric perturbations into independent "even" 
and "odd" components according to 



9nv = 9^(S) + (odd) + (even) , 



(9) 



where g^, v {S) is the usual Schwarzschild metric. In the 
appropriate gauge, Regge and Wheeler found that the 
odd perturbations may be written as 
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FIG. 1. Evolution of the cluster in the hybrid-PN approx- 
imation. The mean radius of the 20,000 test-particles in the 
initially spherical swarm is plotted as a function of the tidal 
companion radius in the top figure and as a function of coor- 
dinate time in the bottom figure. 
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FIG. 2. The adopted coordinate system. The black hole M 
is at the origin and the companion m is on the z axis. 
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Likewise, the even perturbations may be written, after 
simplifying gauge transformations, as 



M „(even) = 

f ff (r)(l-2M/r) 

o 

V 



ffi(r) 

(l-2Af/r) 







r 2 AT(r) 




r 2 sin 2 e_ff(r) / 

F LO (0) X e (^*) . 

(11) 

We now proceed to solve Einstein's equations to first or- 
der in the perturbation functions as in Zerilli , using a 
delta function (point) source for the companion. We first 
find solutions to the homogeneous equations away from 
the companion and then use the inhomogeneous source 
terms to match these solutions. Setting lu = to comply 
with our quasi-static approximation, Zerilli's equations 
reduce to 



hx = Hx = 0; 



H 



B, = B 



(12) 



d 2 h 
dr 2 
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2r{r - 2M) 



(H-K) = 8irA L . 



(15) 
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In the above equations, Ql, Al, and A£' represent com- 
ponents of the stress-energy tensor. For a point source 
companion moving along a Schwarzschild geodesic, these 
terms are given in Appendix E of Zerilli [|l2| . From Zer- 
illi's tabulation of the source terms, one sees that all the 
source terms go to zero in the static limit except for A^ , 
i.e. Too- 

Homogeneous solutions may be found analytically 
for both the even and odd perturbations. Regge and 
Wheeler show that, in the large r limit (where the com- 
panion resides), there is one solution for ho varying like 
r~ L , and another varying like r L+1 . Since there is no 
source term for h (Ql = 0), the only way that both 
ho and its first derivative can be continuous and regular 
everywhere is to require ho = 0. Therefore, there is no 
contribution from odd parity perturbations. 

The static solutions for H and K are found both by 
Regge and Wheeler jll| and by Zerilli [|l2| in terms of 
hypergeometric functions. Wc arc interested in the low- 
est order nontrivial contribution, the quadrupole piece, 
generated by the companion, so we restrict our atten- 
tion to the L = 2 perturbation. Inside the orbit of the 
companion, the solution which is regular at the horizon 
is 



H 
K 



Kir(r - 2M) 
Ki(r 2 - 2M 2 ) 



Kir 

o 

K±r 



(16) 



where the last equality in each equation above holds at 
large r 3> M and ki is a constant to be determined 
by matching at the companion. These equations may 
be verified by direct substitution into Eq. ( p"4| ) or (|T5|). 
Outside the orbit of the companion, the perturbation is 
given by the solution regular at infinity, 



H 



K-2 



K 2 r 



K 



2M(2M 3 +4M 2 r-9Mr 2 +3r 3 ) 
(2M-r)r 

3r 2 (r-2Af) 2 ln(l-2M/r) 
+ (2M-r)r 

-3 



AM 6r 3r 2 ln(l-2M/r) 



(17) 



+61n(l - 2M/r)} 



K 2 r 



where the last equality in each equation again holds at 
large r S> M and k 2 is another constant to be determined 
by matching. 
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Now we determine the two constants by matching the 
solutions at the radius of orbit of the companion, R. Be- 
cause A®\ the source term of (|l4|), does not go to zero 
in the quasi-static case, n\ and k,i do not vanish, as was 
the case for the odd parity solutions. Taking fl>Mas 
required by the quasi-static approximation, we use the 
asymptotic expressions in equations ( |l6| ) and ( |l7| ) for the 
matching. Requiring that K be continuous across R then 
yields our first condition 



K\R = K2R 



-3 



K2 = R K\ 



(18) 



The second condition may be found in either of two ways. 
First, one may integrate Equation (Q) across the source 
to obtain a jump condition relating difference in the first 
derivative of K on either side of r = R to the source 
strength. Alternatively, one may compute k,\ directly 
by taking the Newtonian limit of H and matching it to 
the Newtonian tidal potential. Either way, one obtains 
K\ = Thus, /ioo, for example, is 

h 00 = (1 - 2M/r)(r 2 - 2Mr) Kl Y 20 (9) 

= (1 - 2M/r){r 2 - 2Mr)-^(3 cos 2 6 - 1) . 

(19) 

From (|l9|), one can immediately verify that the metric 
derived above has the correct Newtonian limit — the po- 
tential reduces to the classic Hill potential. This solu- 
tion for the even parity solutions was first discovered by 
Moeckel Q. 

The entire metric may now be written out, defining 
P = P 2 (cos9) = 3 cos 2 61- 1: 



Obtaining the geodesic equations in this spacetime up 

2 

to first order in the tidal expansion parameter -7^5- , we 
find the equations of motion for a test particle near the 
black hole to be 



3Mf 2 



M(r - M) 
r 2 - A/ 2 ~ (r + A/) 3 



rP 



+ (r - M)n 2 



I (r — M) 



-3rr 2 + (r + 2M)£ m z 
+(r - M)(r 2 + ArM + M 2 )fl 2 



r 2(r 2 - M 2 )%0r 



4M-2r , 



>r + cos f sin f 



■j^P[-2(M + 2r)9f] 

(3r 2 + 2rM - 3M 2 )9 2 



to dP 1 
" R 3 d0 2 



+ ( r 2 + 2rM _ M 2) gin 2 Q ^2 + <Z_Mf + ^ 



4M-2r ± ■ 



2 cot (9 0<j>-2jSrP(M + 2r) 



-2^%(r 2 + rM - M 2 



(22) 



where dot denotes differentiation with respect to 
Schwarzschild coordinate time, and f2 2 = 9 + sin 2 6(f) 2 . 
To avoid coordinate singularities, it is convenient to inte- 
grate ( |22] ) in Cartesian coordinates. To facilitate an oth- 
erwise tedious transformation, we first rewrite the equa- 
tions of motion in 3-vector form. We make the identifi- 
cation 



ds 2 = 



l-2f) + ^P{r-2Mf 



(1 



2 An 



T Pr 2 



dt 2 
dr 2 



' R 

r 2 + j^Pr 2 {r 2 - 2M 2 ) 



dfl 2 



(20) 

We note that the perturbed spacetime (^0|) furnishes an- 
other example supporting the conjecture of Bekenstein 
p3| that the horizon area of a near-equilibrium black 
hole is an adiabatic invariant (see Appendix). 

To compare with the hybrid-PN treatment, we trans- 
form to harmonic coordinates by replacing the areal ra- 
dial coordinate, r, with the harmonic radius, = r — M. 
Dropping the subscript, we obtain in harmonic coordi- 
nates 



ds 2 



M 



M 



dt 2 



M 



M 



dr 2 + {r + M) 2 dQ 2 



+ ^P[(r - Mfdt 2 + (r + M) 2 dr 2 



R 3 



+ (r + M) 2 (r 2 + 2rM - M 2 )dfl 2 



(21) 



(r -rcf? sin 2 9-r9 2 )-§ f 

i.2 



+ (2f6-r<f) 2 sin 9 cos 9)/r)-§g (23) 
+ (<j> + (2r4> sin 9 + 2r6j> cos 6)/(r sin 6)) ^ , 

a vector which we have constructed to be the same as 
the Newtonian 3-acceleration in spherical coordinates 
and we know transforms to (x, y\ z) in Cartesian coor- 
dinates. Identifying the velocity 3-vector v = (f, 9, <j>) in 
the spherical coordinate basis and setting r = r-3- allows 



us to write (123) in the compact form 3-vector form 



Bv] 

m 



+ /3v + 7e r 



(24) 



where A and B are again given by Eqs. (g) and (||). The 
quantities a, (3, and 7 are given by 
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a = (r 2 + 2rM - M 2 ) 



1 • 2 

v — r 



.2 , (r-M) 2 



(r + M) 2 



(3 = -2(3 cos 2 6 - 1)(M + 2r)r 



-12 cos 9 sin 6(r 2 + rM - M 2 )Q 



7 = (3cos 2 6- 1) 



2M(f 5 



(r-M) 2 , 
(r + M) 2 ' 



2 • 2 

+A/(r 2 - 2rM - M 2 ) - ~J 



-12rMr6cos0sin6. 



(25) 

where is the angle between r and R. It is now 
straightforward to write out the Cartesian components 
of Eq. (|24j). To incorporate the (slow) inspiral of the 
companion in the context of our quasi-static approxima- 
tion, we treat R = R(t) in Eq. ( pij ) as a parameter, 
which is slowly evolved in accord with Eqs. (0), (R) and 
#■ 

The leading terms independent of the companion m 
are identical in both (|2~i| ) and (Q), a result which is con- 
sistent with the adoption of hybrid equations in the PN 
analysis. The tidal term in (gj) reduces to the lowest- 
order Newtonian expression used in (tth in the weak-field, 
slow- velocity limit where v <C 1, M/r -C 1. Since the 
test particles orbit close to the hole, the relativistic tidal 
expression will indeed cause departures from the motion 
predicted by the Newtonian tidal term. 



B. Numerical Implementation 

A spherical swarm identical to the one evolved with 
the hybrid-PN equations of motion was evolved with 
Eq. (p4[) . Because collapse occurs later in this simula- 
tion, the swarm was evolved until R/M = 26.30. The 
results are summarized in Figures 3-5. In Figure 3, the 
mean cluster radius is computed once again and com- 
pared to the hybrid-PN result. The qualitative nature of 
the evolution is the same, and the existence of a crushing 
instability is evident. The main effect of the new rela- 
tivistic terms is to delay the collapse until the perturber 
gets somewhat closer to the hole. In particular, of 20,000 
particles, 16,324 fell into the hole (81.6 %) by the time the 
companion reached R/M — 26.30, a slightly smaller per- 
centage than in the hybrid-PN simulation. Once again, 
the compact cluster is observed to be stable in isolation, 
but driven to collapse by the presence of a sufficiently 
strong tidal perturbation. 
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FIG. 3. Evolution of the cluster in the Schwarzschild 
perturbation treatment. The mean radius of the 20,000 
test-particles in the initially spherical swarm is plotted as a 
function of the tidal companion radius in the top figure and as 
a function of coordinate time in the bottom figure. In the top 
figure we compare the results found for our two treatments; 
the behavior is qualitatively similar, but the crushing effect 
occurs slightly later in the fully relativistic, Schwarzschild per- 
turbation treatment. 
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FIG. 4. Snapshots of the cluster swarm at selected times during the binary inspiral. In A we show the view looking down 
along the «-axis, which is perpendicular to the orbital plane of the companion; in B we show the view looking along the x-axis, 
which lies in the orbital plane of the companion. For the top frame the companion is at R/M = 60; for the middle frame the 
companion is at R/M = 30.7; for the bottom frame the companion is at R/M — 26.30. 
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As the geodesic equations in the perturbative rela- 
tivistic treatment are derived from a self-consistent La- 
grangian, they satisfy strict conservation laws even in 
the presence of a companion, assuming it is stationary. 
These conservation laws provide a means of testing our 
code and particle integration scheme. For example, if the 
companion is fixed at an arbitrary position on the z-axis, 
the perturbed Schwarzschild spacetime ( ^l| ) admits two 
Killing vectors, Jj and yielding conservation of par- 
ticle energy p t and angular momentum p^ >1 even for large 
tidal fields. Given that we linearize the tidal field and 
retain only the lowest order terms in m in our equations 
of motion (p4|), energy and angular momentum conser- 
vation arc no longer exact. However, we have tested our 
code and have shown that it reliably obeys these conser- 
vation laws to the required order. 

The evolution of the swarm is depicted in Figure 4, 
where snapshots of the cluster are shown at three dif- 
ferent times from two different viewing angles. By the 
end of the simulation, the original spherical swarm is re- 
duced to a sparse cylindrical band of particles whose axis 
is perpendicular to the orbital plane of the companion. 
Apparently, particle orbits at small inclination angles to 
this plane are more stable than those which are perpen- 
dicular, a result already noted in This feature is 
evident in Figure 5, which plots the trajectories of three 
representative particles. The simulation confirms that 
gravitational collapse can be induced in a collisionless 
cluster by the tidal field of a binary companion. 



APPENDIX: ADIABATIC INVARIANCE OF THE 
HORIZON AREA 

Bekenstein |l3| has argued recently that the horizon 
area of a near-equilibrium black hole is an adiabatic in- 
variant. Using the focusing equation, he derives the fol- 
lowing conditions which must be satisfied in order for 
area to remain fixed: 

T af3 ri p = 0; C a/ 3 lS l a mPpm S = , (on the horizon) 

(Al) 

where C a p 1 s is the Weyl conformal tensor, and l a and m a 
are legs of the Newman-Penrose tetrad (l a are tangents to 
the null generators of the event horizon). For our quasi- 
static scenario in which the companion undergoes a very 
slow inward spiral, there is no matter and and esse ntially 
no wave flux at the horizon, so both conditions ( Al) hold. 
(Other situations in which adiabatic invariance has been 
demonstrated are discussed in [H| and fll8|.| We may 
use our metric in Schwarzschild coordinates ( [20| ) to ver- 
ify explicitly that in accord with Bekenstein's conjecture, 
the horizon area remains invariant in our model as the 
companion is brought in slowly from large distance. 

First, rewrite ( J20| ) in Kruskal-Szekeres coordinates, ig- 
noring angular terms, reintroducing areal coordinate r 
and defining k = -gg-, 




FIG. 5. Test-particle trajectories about the black hole 
during the inspiral of the binary companion from R/M — 60 
to R/M = 26.30. Frame (a) shows the orbit of a particle mov- 
ing in the orbital plane of the companion; it is not captured 
by the time the integrations terminate. Frame (b) shows the 
capture of a particle initially orbiting in a plane perpendicu- 
lar to the companion plane. Frame (c) shows the precession 
and capture of a particle initially orbiting at an angle of 45° 
to the companion plane. 
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dg2= 32A£3 e _ r/2M f _ dv 2 +du 2 



+ - 



2Pnr 



-r/2M 



(u 2 + v 2 ){du 2 + dv 2 ) — 4uu du dv 



(A2) 



Note that this metric agrees with that derived by 
Vishveshwara |L9| for general even perturbations of 
Schwarzschild in Kruskal coordinates. Null geodesies are 
found by setting ds 2 = 0, for which 



du 

dv ~ 2(1 + (u 2 



8uvKPr -r/2M 



M 



2 \ 2kPt 
M ' 



V 2 ) 



r/2Af) 



± 



ivnPr) 2 — r/M 



( u 2 + u 2)2(2^)l e - r /Af_ 1 



(A3) 



Direct substitution verifies that (A3) is satisfied by u — 
±v, r(t) = 2M. Therefore, the horizon on the perturbed 
black hole remains at 2M. The surface area of a shell of 
constant radius is given by 



A = I I V9~e 
>o Jo 



(A4) 



Due to the orthogonality of Legendre polynomials (P2 
and Pq in this case), the quadrupole perturbations in 
the angular piece of metric ( |20| ) do not change areas of 
shells. Therefore, the area of the spherical shell r = 2M, 
the location of the horizon, remains constant. 

Hartle |2(| has studied the effect of tidal fields from 
more general stationary sources on the horizon area of 
slowly rotating black holes. In the limit of a nonrotat- 
ing black hole, his analysis gives a vanishing time rate of 
change for the area, consistent with our explicit calcula- 
tion above. 
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